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Abstract 


This report covers progress on NASA Grant NSG-2010 during the period 
September 1, 1975 to February 29, 1976, The coding of our 3-D upper atmo-f 
spheric model for the ILIAC computer at NASA Ames Research Center is now 
nearing completion. Dr. Alyea presented a paper on the 3-D model at the 
Australian Conference on Climate and Climatic change held on December 7- 
11, 1975. Dr. Prinn attended the Space Shuttle Environmental Workshop on 
Stratospheric Effects in Houston, Texas on March 24-25, 1976 and reviewed 
our progress on modelling the ozone perturbations resulting from the Shuttle 
booster exhaust. Mr. Gary Moore, a graduate student, has been studying 
a time-dependent version of our 2-D model and is also investigating the 
sulfur cycle in the stratosphere. Drs Alyea, Cunnold and Prinn have recently 
been emphasizing the significant role of meteorology in influencing strato- 
spheric composition measurements and in this regard Dr. Prinn presented a 
paper on this subject to the NASA Stratospheric Advisory Committee on Dec. 

17, 1975 at NASA Goddard Space Flight Center. This paper will shortly appear 
in the Bulletin of the American Meteorological Society. 
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1. Introduction 


We will discuss significant progress over the six-month period covered 
hy this report under 3 main categories: the 3-D model with full chemistry 

and its coding for the ILIAC IV, the impact of meteorology on measurement 
programs and measurement interpretation, and finally further work utilizing 
our 2“D model. 

2. 3-D Chemical Dynamical Model 

With the help of a full time programmer, the present 3-D model (with 
abbreviated chemistry) is being coded for the ILIAC IV computer and 
this arduous task is now nearing completion. We continue to use our GDI 
Teleterra 1030 terminal to aid this process. In the near future we should 
he ready to run a simple test version of the full dynamical program. The 
chemistry in the ILIAC IV version of our model will contain about 40 chemical 
reactions including those involving formation and removal of chlorine 
nitrate (CIONO^) • Our present full time programmer, M, Batish, will probably 
he replaced sometime this summer with a full time research scientist. We 
should have most of th'e programming complete hy that time and he ready to 
run scientific experiments. We are presently considering suitable candidates 
with a good background in spectral modelling. 

Dr. Alyea presented a paper on some of our 3-D modelling work at the 
"Australian Conference on Climate and Climatic Change" held in Clayton, 
Victoria, Australia on December 7-11, 1975. A copy of the paper entitled 
"Simulation experiments on the climate of the stratosphere and stratospheric 
ozone distributions under natural and man-made conditions" is appended to 
this report (Appendix I). 


Dr* Prinn attended the ’’Space Shuttle Envxromental Workshop on Strato- 
spheric Effects” held at the NASA Johnson Space Center on, March 24-25, 1976. 
He gave a brief review there of our plans for simulating the ozone deple- 
tions resulting from the booster exhaust products. Some prelimimary runs 
of our 2-D model have been made and relevant results will appear in the 
workshop report. We do not have any de"inite numbers to report at this 
time and indeed we feel our best estimates will be made only after we have 
in operation the ILIAC version of our 3-D model with full chemistry. 


3. Stratospheric Dynamics and Stratospheric Measurements 


We have been particularly concerned over the past year with the problem 
posed by dynamically induced stratospheric variability and its influence on 
the planning and interpretation of measurements for minor ccastituents. 

Dr. Prinn presented a paper on this subject ot the NASA Stratospheric Advi- 
sory Committee Meeting held at the NASA Goddard Space Flight , Center on. 12/17/73. 
This paper entitled "The impact of stratospheric variability on measurement 
programs for minor constituents" is appended to this report (Appendix II). 

It has been accepted^ for publication in the Bulletin of the American Meteoro- 
logical Society and will appear in the June, 1976 issue. We feel that a 
correct interpretation of stratospheric minor constituent measurements will 
entail careful use of 3-D modelling results and the kno^<m meteorological 
data for the observing sites. We also feel that observers should coordinate 
their observations in a much more meaningful manner than presently ex'sts. 

We arc presently studying a fex-j examples of recent fluorocarbon and tCl 
measurements in order to illustrate these ideas. 
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4. 2-D Chemical Model 

Mr. Gary Moore, a graduate student, ties recently. completed a study of 
the effect of continuous time updating of. transport parameters in 2-D models. 
He has written a detailed report entitled ”A partial report on the time 
updating modifications made on the 2-D Stratospheric model, 2AM4" which will 
provide us with a useful reference for ftiture use of this modelling approach. 
Mr. Moore has also heen studying the cycle of sulfur in th& atmosphere with 
particular emphasis on the production of stratospheric H^SO^ from the GO S 
and SO^ derived from fossil fuel comhustion. He intends to ultimately model 
the stratospheric sulfuric acid layer utilizing our 2-D model hut adding 
in the sulfur chemistry and modelling the important tropospheric removal 
mechanisms (rainout, washout, etc.) in some detail. 
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ABSTRACT 

The scales of temporal and spatial variations in stratospheric minor 
species concentrations and dynamical quantities are "briefly reviewed. De- 
spite the fact that the stratosphere is often considered as a very quies- 
cent region of the atmosphere » it exhibits transient and generally unpre- 
dictable episodes of intense activity. The resultant necessity for some 
degree of coordination between different observers in their choice of 
locations and times for measurements of minor constituents is emphasized. 
¥e provide the specific suggestion of one or more short "International 
Stratospheric Periods" whose duration will depend upon the time of year. 
During these short periods as many observers as possible would perform, an 
integrated series of observations designed to optimize the usefulness ■ of 
their data to atmospheric scientists. 
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1 , Introduction 
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Over the past few years, atiaospheric scientists have hecome increas- 
ingly aware of the necessity and urgency -for measurements of the concentra- 
tions of a significant number of stratospheric minor constituents. Of 
particular concern has been the need for a sound observational basis for 
evaluating the stratospheric ozone models that have recently been used for 


predictions of ozone depletion by anthiopogenic trace gases such as NO and 
NO 2 from supersonic aircraf tj and CFCl^ and CF^Cl^ from aerosol cans and 
refrigerators. Considerable effort has been taken to delineate those spe- 
cies and dynamical variables whose measurement will contribute most to these 
model validations. In contrast, the possible need for a significant degree 


of coordination between different observers in their choice of times and 


locations for these measurements has received comparatively little attention. 
For species with chemical time constants which appear to he smaller 


than typical advection times, there are a number of persuasive arguments 
for considering simultaneous measurements as discussed at a recent workshop : 
on the stratosphere (NASA, 1975). A few of such types of observations in 
fact already exist (Ackerman et al, 1975; Lowenstein and Savage^ 1975; Evans 
et al, 1975). Examples of groups of species for which siraultaneous measure- 
ments are highly desirable are: odd nitrogen, (NO, NO^, HNO^) , odd oxygen 

(0, 0^^’ chlorine (Gl, CIO, HCl), and odd hydrogen (OH, HO^, ^2^2^* 
many of these short-lived species are predicted to exhibit significant 
variations with time of day (Whitten and Turco, 1974) and with latitude 
and season (Frinn et al, 1975), the hour and location of such observations 
is clearly an important ingredient in their interpretation. Fortunately 


J 


most observers appear fully aware of such Goiisideration,s so we need not 
. pursue these requirements any further here. 

However, somewhat less emphasis has heen given to the coordination 
of Tneasurenients involving species whose concehtf at ions depend pn Both chem— 
ical and _advective processes . ^Utilizing observations and .numerical model 
results we can readily show' that the temporal and spatial variability of 
such minor species is significant. In addition, we know there exist impor- 
tant variations in stratospheric temperature Cwhich will influence the 
chemical terms) and stratospheric motions (which will govern the advective 
terms). Such inhomogeneities should obviously be fully appreciated rn\ 
planning o'^servational programs. ’ * 


Hynamical variability 



let us begin by looking at a few examples' of the pseillati OHS in 
namical variables in the stratosphere. In doing so we will- tend to choose 
cases wliich will -illustrate the*~greatest~variahi Should -not 

. be-interpreted-as an attempt -to over-state our -nase but -rather as pro-- 


viding a review of the significant changes which can: occasionally occur 
in the stratosphere. 


Rapid alterations in stratospheric temperature are of course well 
knot™ to meteorologists. The most Impressiye are the Sudden Stratospheric 
Warmings and in Fig. 1 we show the temperature changes observed over Scot- 
land (57 “N) during the 1967-68 warming (Quiroz, 1971) .: Of particular in- ; 
terest to our discussion here is the fact that, the temperature increased 


by :8Q°C over a ten-day period at the 43 km. ^1^^ This pa.rticiilar vjarming, 

was indeed strong enough to cause the usual winter circumpolar vortex 
(westerlies) to be replaced temporarily with anticyclonic easterlies. Such 
impressive events appear to occur about once every one or two years. 



■ ’ .■ -4-. ■ ■■ 

Warmings on a smaller scale octiur with greater ftequency and ill a gen- 

erailly unpredlctaBle manner. The magnitude and frequency- of these • . 

. minor .warmings is particularly well illustrated in Fig*. 2 in which . 

teftperatures at the 30 mb level at high and mid— latitudes are given for 

two winter seasons (Madden, 1975)1 The observed temperature oscillations 

typically exhibit 1-3 week periods. The 1967-68 major warming is also 

clearly illustrated. In recent unpublished runs of oiir 3‘-dimensipnai 

stratospheric circulation model at M. I, T. (see Cunnold et al, 1975), we 

have also obtained many minor winter warming events and these are shown 

as functions of longitude and time at 45 km in. Fig. 3. In addition to 

the general warming at 45 km the reader should note what appears to be 

' -1 

a Rossby wave with zonal wave number 1 travelling westward at 3 ra sec 
This produces temperature oscillations ^^ 50“K in addition to the general 
warming. At 30 km specific wave numbers are less obvious hut significant 
temperature oscillations are present there also. As a further example, 
northern hemisphere temperature changes observed over a 7-day period in 
March, 1964 are shown in Fig* 4 (Efewell et al, 1974). Alterations in 
vertical velocities deduced: using the thermodynamic equation are also 
depicted in the: latter illustration. 

The standard deviations for the stratospheric meridional wind during 

, - , - ■ -1 

the winter imply typical variations ranging up to 25 m sec . Standard 
deviations during the summer are. usually about a factor of 3 or 4 less 
(Mewell et al, 1966). The daily variation of the zonal wind during March, 

1964 is show in Fig, 5 (Dopplick, 1971) . At 30 km altitude, 15 m sec ^ 

■ 

westerly -td-nds reverse to 15 m sec easterlies in a period of ' about 2 
Weeks. For - estim^ daily variations in advective terras we are of course 
cdncerned with daily correlations between the wind field components and 



the various transported quantities. The variahility of the northward 
eddy flux of heat, at 60“K can he equated with the suhstantlal oscillations 
in the convergence, of this heat flux northward of 60°N which are shown in 
Fig. 2. As a further example, the daily variation of the vertical eddy 
flux of geopotentiajL is given in Fig. 6 for February, 1964 (Dopplick, 

1971), The substantial variations in this energy flux with both latitude 
and time should be particularly noted. We will give further examples of 

^ the daily .variability of. suc h cor relal^ons shorty. __ , 

Transport processes across the tropo.pause are also of particular 
concern- for quantitative studies of upper atmospheric trace constituents. 
Reiter (1975) estimates that during one year about 40% of the total strato- 
spheric air mass is exchanged with the troposphere through the tropical . 
Hadley cell, 20% through large— scale eddy transports in the jet stream 
region, and 10% through seasonal changes in tropopause .heights'. . The up- 
ward flux in the tropical branch of the Hadley cell mainly occurs through 
localized and transient cumulus convection (Riehl and Malkus, 1958),; while 
exchange of air in the jet stream region can he related to tropospheric 
cyclonic activity (Mahiman, 1965). Thus these two processes are expected 
to have significant temporal and spatial variations. In Fig, 7 we illus- 
trate. the surprisingly large variahility in the January, 1973: tropopause 
position (Reiter, 1975). The July, 1973 variahility is somewhat smaller 
but nevertheless still significant. This has implications both for . 
tropospheric-stratospheric exchange and for variations in the dynamical 
stability (eg. Richardson number) in the 100-500 mb region. 

The differences between January and July just noted remind us we should 
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also discuss scales of 'variatn.lity during different seasons. Both ttieofy .. . . 
and observation (see Holton, 1975 for a review) suggest tbat yertically- 
propagating extra-tropical planetary waves provide the major energy source 
for driving the lower stratosphere. , This upward flux is modulated in part 
by downward reflection of quasi— stationary Rossby waves by east erly winds 
in summer and strong westerly winds in winter (Cbarney and Drazin, 1961) . 

This modulation is evident in Rig. 8 where we have shown, the daily verti- 
cal eddy flux of geopotential through the 100 and 10 mb levels during 1964 
(D opp lick,: 1971) . Similar annual cycles, in the variance of temperature 
.and meridional wind are evident at most stratospheric levels in unpublished 
data here at M. I. T. We therefore expect the most rapid and largest varia- 
tions in dynamical variables during the. late winter months and indeed most 
of the examples we have just provided are derived from this period. Of ^ 

■particular importance for our discussion here, : is .the time-scale over 
which we may expect no signif icant dynamical changes and this is clearly 
■ very different in summer and winter seasons . 

3. Minor constituent variability 

^ Having thus noted , the dynamical variahilty of stratospheric air , we 
can move' on to discuss si.milar effects: for,, stratospheric minor cons titt-- . ; 
ents whose local concehtrations are affected by the scales and velocities 
of mixing processes as well as chemis try,. Although the ma j or ity of the 
Imp of tant longer-lived stratospheric species have .not been measured at 
. ^ a wide range of. Ibc^itids and :^over 'extended; periods, of timey .enough- data -y 

does exist (eg. for 

the types of iuhpmogeneities we should expect in stratospheric air. .As / 
we implied earlier, our principal thesis will be that uncoordinated 



observations of different species at different localities and/or times are 
in general incompatible. As a result erroneous estimates of vertical 
profiles may result from simple superposition of such uncoordinated measure- 
ments . 

Let us look at some typical early spring ozone variations. In Fig. 9" *. 
we illustrate the temporal and latitudinal variations in total ozone ob- 
served during Marctij 1964 (Newell, et al, 1973) . It is clear that total 
ozone increases of up to 50 % may occur at mid- latitudes over a time span 
as short as ten days. The quasi-horizontal eddy activity sho^ra. in Fig. 9 
is also increasing very significantly during, this particular period. Ob- 
servations at the same latitude but 10 days apart would therefore provide 
significantly different results. Similar rapid temporal variations occur 
in ozone concentrations at various altitudes as illustrated over Switzer-: 
land in Fig. 10 (Dutsch, 1974) .' In the latter case we see that factor of 
B changes over a period of 33 days are not uhusual. Sighif leant variations 
in total ozone (;$ 25%) • also occur with longitude (Reiter, and Lovillj 1974)' 
as illustrated in Fig. 11. for 20°N latitude, Host readers will be familiar 
with the significant variations in monthly mean total ozone as a function 
of position and season (Dutsch, 1971) and such variations have been simu- 
lated in the K.I.T. model (Alyea et al, 1975) . 

Unfdrtuhately, measurements of water vapor in the stratosphere are 
few and generally restricted to Europe and 'North America. Spatial and 
temporal variabilities are certainly expected because of temperature changes 
and the localized and transient injections of water vapor into the lower 
stratosphere associated with the tropical up.welling (ie.> cumulus towers) 
and with mid-latitude thunderstorms. , We can illustrate our case using . 
recent observations of the latitudinal dependence of the H^O mixing ratio 


in the 15 - 20 Ion region during October, 1973 through January, 1974 whic*. 
are illustrated in Fig. IZ (Kuhn et al, 1975). Note the doubling of this 
ratio in the region of the ITCZ. 

Observations of hydrogen nitrate and sulfuric acid have exhibited 
significant altitude and latitude variations, A stmmary of global observa- 
tions during Northern Hemisphere Spring by Lazrus and Gandrud (1974) is 

given in Fig, 13 and the HNO^ observations can be compared with our 2-D 
model results (Frinnet al, 1975) which show qualitative agreement with 

these data. The map is constructed hy superimposing various obser- 

vations during the Spring of 1971 and agreement with a 2-D model in which 
short-term variations are essentially averaged out is perhaps not surpris- 
ing, The process of superposition is of course contradictory to the gen- 
eral theme of this paper hut our purpose in presenting such results is 
to qualitatively illustrate the spatial variability of HNO^, Significant 
spatial variations suggest that significant temporal variations will 
result during sporadic rapid large-scale interchange of air parcels of 
the type discussed earlier. The fact that lower stratospheric mixing ap- 
pears to he dominated by large-scale rather than very small scale motions 
(Gunnold et al, 1975; Heck and Panofsy, 1975) reinforces this notion. 

As a final example we show in Fig. *14- the observed latitudinal varia- 
tions in CClgF during April and October, 1974 (Krey and Lagomarsino, 1975), 
Variations are most significant in the vertical and temporal variations over 
at least a 6 month period cannot be neglected. Similar results were obtained 
by these workers for SF,. 


4. An International Stratospheric Period 

It is hoped that we have provided a suff.icient case for concern over 
the inhomogeneity of the stratosphere. . We particularly emphasize that 
superposition of profiles of different species measured at different posi- 
tions and times is in general hazardous. In addition, we must he prepared 
to accurately scale localized vertical profile ohservations to the annual 
global average for comparison ^tri-th 1-D models in which most of the ozone 
predictions have been made. Multi-dimensional models, particularly 3-D 
models where dynamics are predicted, will play an important role here hut 
we are still faced with the problem of day to day predictability of the 
stratosphere. We must therefore obtain some idea of temporal variations 
on a monthly or seasonal scale. These observations need not be continuous. 
i.n space and time providing imaginative use of 3-D numerical models and 
existing stratospheric observations is made. 

With these points in mind we would like to recommend the designa- 
tion of one or more "International Stratospheric Periods" (ISP) for the 
purpose of coordinating existing measurement projects. In each ISP we 
suggest the following scenario for the timing and location of observations 
and their interpretation: 

1) Each ISP should be defined with enough lead time being provided 
to observers to prepare for participation,. The existing funding agencies 
in North America, and if possible in Europe and Australasia, should en- 
courage their various observing groups to cooperate fully during these 
periods, 

2) The duration of each ISP will vary with the time of year chosen. 

♦ 

Based on our earlier discussion we suggest intervals of about 1-3 days 



during late Fall, Winter, and early Spring; about 1 weekf during early 
Fall and late Spring; and about Z weeks during Summer. If one ISP is 
planned each year then a different season could be chosen each time, 

3) Balloon and rocket flights during ISP should be as far as pos- 
sible chosen so as to ensure measurements of vertical profiles of as 
many constituents and d 3 maraical variables as possible above each station. 

In the event that this is not possible, superposition of the results from 
different stations at the same latitude and local time of day is generally 
preferable to superposition of results from stations of the same longitude. 
The aim would be to obtain self-consistent sets of vertical profiles for 
at least one high-latitude, one mid-latitude and one low-latitude station. 
Southern hemisphere coverage would initially be sparse but would hopefully 
improve in later years, 

4) Aircraft flights and satellite observations during ISP should 
cover the regions above and between the stations in 3). This will pro-' 
vide knowledge of horizontal variations and help to relate the localized 
measurements at each station to the global picture at the time of ISP. 

Some coverage of the Southern hemisphere would also be useful- for assess- 
ing the global picture. 

5) As soon as possible after the ISP all the results should be com- 
piled Into a compact data report made readily available to the atmospheric 
scientific community. This report would be recognized as preliminary; it 
would not preclude subse<iuent publication of both the data and methods in 
the open literature. Recognizing the forcing role of the troposphere, this 

volume should also include pertinent global tropospheric and lower strato- 
spheric meteorological data obtained during the ISP under existing daily 
atmospheric observation programs. 


5. Concluding remarks 

¥e feel tiiat a scenario of the type proposed here will considerahly 
increase the usefulness of future stratospheric measurements for aeronoraers 
and meteorologists alike. To conclude, we should emphasize that except 
for funding of a small staff for coordination and data report compilation 
purposes, there is no requirement for additional funding over and above 
that provided by existing agencies. There should be no hint of a central 
controlling agency; each observer would have to obtain his oxto funds based 
on the 'merits of his particular experiment. T^e also recognize that the 
vagaries of instrument development and the' weather may prevent participation 
by some observers during a particular ISP. Hox<rever, even partial success 
in such a venture must he deemed superior to allowing a haphazard and 
totally uncoordinated program for measurement of stratospheric minor con- 
stituents to exist. 
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Figure Captions 


1. Stratospheric temperatures T (°C) observed over West Geirinish, Scot- 
land (57°N) at various altitudes (km) dur 7 .ng and after the December 3 1967 
Sudden Stratospheric Warming (after QuiTroz, 1971). 

2. Zonal mean temperatures T (°C) for 70-80°KT (HIGH LATITUDES) and 
40-50°N (MIDDLE LATITUDES) during the winters of 1967-1968 (upper) and 
1971-1972 (lower). Also illustrated is the convergence of the meridional 
eddy heat flux (-G/day) northward of 60°H (EDDY COHV) . The latter flux 
is defined as the zonally- aver aged correlation 'between the geostrophic 
perturbation meridional velocity v* (m sec”^) and the perturbation tem- 
perature T* (“C) j with only the first six zonal wave numbers being con- 
sidered. All data refer to the 30 mb level. 

3. Stratospheric temperatures T (°C) obtained during January and February 
at 70“N from a run of the MIT Stratospheric Circulation Model. Data are 
given for the 2 mb level, 

4; Northern hemisphere stratospheric temperature T (°C) observed at the 

10 mb level on March 19 and 12 during the 1964 spring warming (solid lines) . 

—5 —1 

Also sho’Eim (dashed lines) are vertical velocities u (10 mb sec ) computed 
from observations for the same period using the thermodynamic equation (after 
Newell et al, 1974). 

5. Daily longitudinally-averaged zonal winds u (m sec ) at the 10 mb level 
computed from Northern hemisphere observations assuming geostrophy for March, 
1964 (after Dopplick, 1971) . 


6. Daily zonally-averaged correlations in the Northern Hemisphere between 

“5 ”1 

vertical perturbation velocity in the Northern Hemisphere W* (10 mb sec ) 

and perturbation geopotential height Z* (cm) at the 10 mb level during 

February, 1964 (after Dopplick, 1971). This correlation has the units of 

“2 ' 

energy flux (deciwatt ,m ) and constitutes the pressure work contri- 

bution to the total vertical wave energy flux. A negative correlation 
denotes an upward flux, 

7« Cumulative percent frequency distribution of tropopause pressures 
over North America during January, 1963 . (after Reiter, 1975). 

8. Daily horizontally-averaged correlations in the Northern hemisphere' 

“3 -1 

between w'* (10 mb sec ) and Z* (cm) at the 10 and 100 mb levels during 

-2 

1964 (after Dopplick, 1971). The units are deciwatt ra (see Fig. 5 
for physical interpretation) . 

9. Total columnar ozone (cm, mb) and’ zonally averaged correlation 
between w* (mb sec ^) and v* (m sec ^) at the 100 mb level during 1964 
Northern hemisphere March, The latter correlation is positive 'for north- 
ward and downward (ie. quasi-horizontal) eddy transport (after Newell, et 
al, 1973). 
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10. Ozone partial pressures (10 mb) and atmospheric temperatures (”C) 
observed as functions of pressure (mb) over Payerne, Switzerland on Feb- 
ruary 20 (solid lines) and 23 (dashed lines), 1970. Lines of constant 
ozone mass mixing ratio (10 gm 0^ per gm air) are also show (after 
Dutsch, 1974). 


11. Total columnar ozone Ccm. mb) observed at 20*’N during June 

1969 by the NIMBUS III IRIS experiment (after Reiter and. Lovill, 1974). 

12. Average water vapor mixing ratios (10 ^ gra H 2 O per gm air) observed 
in the 15— 20 Inn altitude region during the period October, 1973 through 
January, 1974 (after Kuhn et al, 1975). 

13. Aircraft 'observations at various altitudes (km) of the mass mixing 

—9 

ratios (10 gm per gm air) of HNO^ gas (lower graph) and sulfuric 
acid aerosol (upper graph) during the 1971 Northern hemisphere Spring 
and Southern hemisnhere Fall (after Lazrus and Gandrud, 1974). 
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14. Volume mixing ratios (part per 10 ) of the chlprofluorocarbon 

CCl^F observed at various altitudes (km) during April (left hand side) 
and October (right hand side), 1974 (after Krey and Lagomarsino, 1975). 
The tropopause position is indicated by a dashed line. 
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Although ozone (0^) constitutes less than one-millionth of the mass 
of the atmosphere, it plays an important role in heating the stratosphere 
and in controling the amount of ultra-violet (UV) radiation that reaches 
the ground. The global distribution of 0^, however, evolves from a com- 
plex mixture of dynamical and chemical processes which lead to considerable 
spatial and temporal variations. Thus, in order to simulate the climato- 
logical cha^jges in the stratosphere, a combined dynamical chemical three- 
dimensional model, global in scale, is desirable. 

We are presently using such a model, 1) to simulate the large-scale 
climatological features and ozone distributions of the stratosphere, 2) to 
estimate the possible effects of anthropogenic pollution sources on the 
stratosphere and UV radiation levels at the surface, and 3) to study the 
dynamical and transport mechanisms relating to the stratosphere. To date 
some progress has been made relative to the first two of these points. 

Work has just begun on the third. 

Model description 

In order to accomplish these tasks, our strategy was to develop a 
fairly simple dynamical stratospheric circulation model (SCM) coupled with 
a relatively few chemical reactions considered most important to strato- 
spheric 0^. The basic dynamical set of equations makes use of a form of 
quasi-geos trophic approximation. The chemistry is limited to the Chapman 


oxygen reactions 


0^ + tiv -> 20 
0 + 0^ + M 0^ + M 
Og + hv 0^ + 0 

0 + Og “> 20^ 

coupled with three odd nitrogen reactions 

KO + 0^ NO^ + 0^ 

NOg' + 0 KO + 0^ 

Kp^ + hV ->■ HO + 0 

For details of the model see Cunnold, (1975). 

Although the model encorporates 26 vertical levels from the surface 
(including orographic effects) to approximately 72km height (wall into the 
mesosphere) , we are not concerned with reproducing the detailed observed 
motions i:i either the troposphere or the mesosphere. The troposphere has 
been included solely to provide a more or less realistic energy source to 
drive the stratosphere, while a number of model levels in the mesosphere 
have been retained in order to remove the upper boundary as far as possible 
from the stratospheric levels of interest. The 3-dimensional prediction 
variables are limited to vorticity, temperature, and ozone although, as in 
many other quasi-geostrophic models, horizontal mean temperatures are not 
predicted but are externally specified. Coupling between the photochemistry 
and the dynamics is achieved not only through the temperature-dependent 
photochemical reaction rates but also through the heating term in which the 
absorption of solar radiation by the predicted ozone is integrated explicitly 
using a simplified radiative transfer law. In the middle and upper layers 




of the model the radiative cooling is approximated using Newtonian rate 
coefficients computed by Dickinson (1973), Below about 20km (lower strato- 
sphere and troposphere) heating is specified in a quite different manner. 

In earlier versions of the model this was accomplished through a simple 
linear parameterization of the Manabe and S trickier (1964) type in which 
the heating rate is proportional to the difference betovfeen an "equilibrium” 
temperature, T*, and the local model predicted temperature. Spatial and 
temporal distributions of vjere determined from data obtained from the 
M.I.T. General Circulation Project. Currently, however, we are replacing 
this linear heating mechanism with total heating data obtained from the 
results of a one year integration of the 13-layer general circulation model 
(GCM) of the Geophysical Pluid Dynamics Laboratory (GPDL) in Princeton, N.J. 
This procedure has the advantage of providing the model with a fairly de- 
tailed heating input, global in extent, and over the entire annual cycle. 

In contrast, the present heating parameterization is based on fairly smooth 
data which is limited to the Northern Hemisphere. 

Horizontal representation of the model variables is accomplished by 
application of truncated series of spherical harmonics. To date all of 
the model computations have been limited to series representations which 
have been truncated after planetary wave number 6 (with equal resolution 
in the north-south direction) , Higher resolutions are undoubtedly desirable 
but we appeal to the tendency for larger scales in the stratosphere. In 
addition, linear theory suggests that only the very long planetary scale 
waves o^ the troposphere penetrate far enough into the stratosphere to act 
as a major energy source for stratospheric motion. 



Because of a lack of good tropospheric heating data in the Southern 
Hemisphere, it was convenient for the first few model integrations to 
assume that the two hemispheres are geographically similar. Thus, the oro- 
graphic heights were defined hy imagining that the Southern Hemisphere is 
a mirror image of the Northern Hemisphere (Fig. 1) . Tropospheric heating 
values from the Northern Hemisphere are also used in the Southern Hemisphere, 
but with a sir month time lag. In this way we have retained realistic inter- 
actions between the hemispheres while making use of the limited tropospheric 
heating data available. 

Simulation of the unperturbed stratosphere 

The initial model integrations were concerned with obtaining a fairly 
realistic climate simulation of the observed stratosphere. The model was 
integrated for a period of 3 years, starting from an initial state of rest, 
until the 0^ and dynamical quantities arrived at equilibrium solution 
conditions. Figure 2 shows a winter-summer seasonally averaged teraperatur® 
cross-section obtained from year three (season 10) of the results of one 
of these model runs (Run,. 12.). The upper part of this figure contains the 
corresponding obseirvational analysis of Newell (1969) . Throughout the strato- 
sphere the model temperatures appear adequate, with some exceptions near 
the poles. Near the winter pole the model predicts temperatures that are 
approximately 20 degrees too low at upper and middle stratospheric levels. 

At the lowest stratospheric levels ('^10-20km) both polar regions of the 
model exhibit temperatures somewhat higher than observations (by approximately 
10-20 degrees). However, the probability of an error in the observational 
analysis exists at about 50km in the summer hemisphere (Netirell, personal 
communication). Here the observational cross-section shows a closed 273“K 


contour which is about 20 degrees cooler than other independent observa- 
tions and also as predicted by the model. In the lower mesosphere, the 
model fails to predict the observed reversal of the pole to pole tempera- 
ture gradient. This probably results from the artificially imposed upper 
boundary condition in which the vertical pressure velocity is assumed to 
be zero. 

Corresponding cross-sections for the seasonally averaged zonal wind 
field (winter-summer) are contained in Fig, 3, The data is as expected 
from the temperature profiles and the thermal wind relationship. The 
climatological winter- summer mean zonal winds obtained from the model 
exhibit the essential lai'ge-scale characteristics of stratospheric motion. 

A winter-hemisphere upper tropospheric (lower stratospheric) westerly jet, 
with maximum Tvind speeds of 40 m/s, is found at an altitude of 10km. in the 
subtropics (although about 15“ too far toward the equator) , In the upper 
part of the stratosphere, a westerly T-jind belt (polar night jet) exists 
in winter, while easterlies prevail during the suirmier. The failure for 
the model upper stratospheric wind belt contours to close with the height 
above the '^601cm level is again attributed to the artificial upper boundary 
condition. 

Another important climate parameter, and one that is certainly impor- 
tant to the transport of stratospheric ozone, is represented by the meridional 
mass circulation. Calculations of this quantity from atmospheric wind data 
have been made by Louis (1974). Figure 4 shows Louis’ mass meridional 
circulation pattern for the winter-summer months, December through February. 

We see that the troposphere contains six individual circulation cells domi- 
nated by the direct (Hadley) tropical cell of the winter hemisphere. How- 
ever, above the tropopause there is a rather abrupt change to a three cell 
(possibly four cell) pattern. Here again, the direct circulation jattera 
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assocxated ’c<rith the tropospheric winter hemisphere Hadley cell dominates, 
but has now eKtended its influence throughout the subtropics and into the 
midlatitudes of both hemispheres. The corresponding mass circulation pattern 
for the months December through February from the model calculations are 
contained in Figure 5. The results are essentially as computed by Louis, but 
the direct tropical cell in the summer hemisphere extends upwards through 
the entire stratosphere instead of being confined to the loiter stratospheric 
levels noted in the observations. Nevertheless, the model exhibits a large 
stratospheric mass transport across the equator from the summer hemisphere 
to the winter hemisphere. This is in accord with the observational evidence 
and thus provides a satisfactory basic meridional flow field for the initial 
poleward transport of stratospheric 0^ away from the equatorial production 
regions. 

Model ozone distributions 

Coupled xflth the dynamical aspects treated above, the model also pre- 
dicts three-dimensional 0^ distributions. In addition to transport of 0^ 
by the wind fields, the model 0^ quantities are also effected by temperature 
dependent chemical reactions. Feedback to the dynamics occurs through the 
heating terms which depends upon the absorption of solar radiation by the 



An important contribution to the quantity of 0_ which exists in the 
stratosphere is realized through the presence of odd nitrogen compounds. 

NO reacts with ozone to destroy it by producing NO^ and O 2 , but the NO^ is 
quickly converted to NO again. Thus, a catalytic reaction chain is realized 
in which the net result is 0^ depletion while the amount of NO remains un- 
altered. This add nitrogen reaction chain has only recently been suggested 
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as being important to stratospheric 0^ distributions and its absence from 
earlier model calculations undoubtedly accounts for their overprediction 
(by about a factor of 2) of 0^. 

Since our model does not -explicitly predict odd nitrogen, we are forced 
to provide spatial and temporal odd nitrogen distributions from independent 
sources. For the model results to be discussed here, we have adopted NO^ 
values as predicted by a two-dimensional model calculation of Hesstvedt 
(1974), from which NO is obtained diagnostically. This NO^ distribution, 
at its midsummer and midwinter extremes, is sho™ in Fig. 6. The values 
lie vrithin the range of uncertainty (about a factor of 2) of the limited 
number of measurements, and we therefore consider them adequate as a first 
estimate. We note from the figure that the NO^ mixing ratios appear to be 
relatively independent of latitude, but vary considerably with height, at- 
taining maximum values at about 35km. 

Figure 7 includes mean zonal 0^ cross sections for the winter-summer 
(December through February) season from observations (Wu, 1973) and from 
the model calculations. The model correctly predicts the do^raward and 
poleward transport of ozone, with maximums at polar latitudes. However, 
in the model, these maximum densities are located about two kilometers 
too high. This undoubtedly results from an inadequate downward transport 
of 0^ into the troposphere by eddy motions which are restrained by the 
modei^s limited resolution (truncated after planetary wave 6). As another 
basis for comparison of observations with our model 0^ results, we refer 
to Fig. 8 in which total columnar 0^ values, as a function of latitude 
and season for the Northern Hemisphere, are shox-m. An observational compila- 
tion by Dutsch (1971) is presented in the upper part while a similar picture. 
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averaged over year 3 of the model results, is contained in the lower pairt 
of this figure. The model has correctly simulated the observed polar spring 
maximum as well as the relative minimum at the pole in the fall. The polar 
spring maximum, however, is some 10% less than is observed, reaching just 
over 400 Dobson units compared to more than 440 in the observations. 

Ozone reduction, by SST^s 

Since the amount of 0^ in the stratosphere is somewhat controled hy the 
presence of odd nitrogen compounds, it stands to reason that aircraft, whose 
engines produce odd nitrogen compounds, flying in the stratosphere may inject 
quantities of odd nitrogen sufficient to be of importance to the 0^ balance. 
This is of particular concern because of the relatively long time periods 
required for the net removal of odd nitrogen at these altitudes. Such long 
residence times are primarily due to the inherent vertical thermal stability 
of the stratosphere. The principle removal mechanism for odd nitrogen in 
the stratosphere appears to involve its conversion to nitric acid (HNO^) ^ 
with subsequent transport into the troposphere and rainout. Thus, it is 
important that we obtain an estimate of the long-range effect on the 0^ layer 
that might be generated by a large fleet of supersonic transports (SST's), 
before such a fleet becomes a reality. 

To do this, we have introduced an additional amount of NO 2 into the 
model to represent a hypothetical SST fleet (Alyea, 1975). This per- 

turbation in was computed using a two-dimensional model for chemical 
stratospheric constituents developed by us at M.I.T. (Prinn, et al, 1975). 

For this first perturbed NO^ calculation, we assumed an odd nitrogen emission 
rate of 1.8 x 10^ metric tons per year, to be released in a corridor centered 



at 45° N and an altitude of 20 km. The corridor extends 15° in latitude 
and 1 km in height. This injection rate was chosen to most nearly repre- 
sent the effect of a large fleet (<^500 aircraft) of the now cancelled 
American SST (Boeing 2707) wh,ich might have existed near the year 2000. 

The results of the two-dimensional calculation of the NO^ perturbation, 
after equilibrium solutions had been obtained (some 25 years of integration 
time), are shown in Pig. 9 (Northern Hemisphere winter). We see that while 
much of the introduced odd nitrogen is found in the Northern Hemisphere, 
interhemispheric transport processes have distributed a significant quantity 
to the Southern Hemisphere, particularly in the 30-40km altitude range. 

The effect of this increased amount of odd nitrogen introduced into 
the stratosphere can be seen in the 0^ distribution simulated by the three- 
dimensional model. Figure 10 shows the seasonally averaged distributions 
of total columnar 0^ as a function of latitude, calculated for both the SST- 
perturbed stratosphere and the unperturbed stratosphere. On a global scale 
the perturbed stratospheric 0^ depletion amounts to about 12%, However, 
the model also shows that about 16% of the Northern Hemisphere 0^ will he 
destroyed. Despite the absence of any direct odd nitrogen injection in the 
Southern Hemisphere, the calculations Indicate that approximately 8% of the 
total 0^ will be destroyed there. The mid-latitude injection of odd nitro- 
gen in the Northern Hemisphere apparently produces a filtering effect on 
the northward transport of 0^ through this region so that there is a larger 
(by about a factor of 2) reduction of 0^ at high northern latitudes than 
at low latitudes. 

To estimate UV radiation reaching the ground, we submitted the total 
0^ distributions for the January and July months to Dr. S. V. Venkateswaran, 
who used this data in a multiple scattering program (N. Sundararaman, ^ al. 



"XU 


1975) to compute daily mean surface UV flux values. The values for both 
the unperturbed and perturbed stratosphere at four wavelengths under clear 
sky conditions were then used to prepare curves of approximate erythermal 
dosage (sunburn producing effect on Caucasian skin). If F(X) represents 
the solar UV flux per unit wavelength arriving at the ground, then the ery- 
thermal dose D. at a given latitude and time of year can be represented by 

U =1 E(X) F(A) dCX) 

A 

where E(X) denotes an erythermal efficiency function (e.g., Cutchis, 1974, 

Fig. 8). T^e results of the numerical integration of this equation for 
the July and January months for which UV fluxes were obtained are shown in 
Fig. 11. The greatest differences between the calculations for the per- 
turbed and unperturbed states occur during July in the Northern Hemisphere, 

The erythermal dose at 30“N for the perturbed state is about 30% greater 
than that for the unperturbed state. At 60“N this difference increases to 
about 58%. Thus, during July the calculations indicate that Northern Hemi- 
sphere mid-latitude erythermal doses under the perturbed conditions are 
roughly equivalent to unperturbed levels some 15° of latitude to the south. 

At the same time, Southern Hemisphere July mid-latitude erythermal radia- 
tion levels are relatively uneff acted. During January, however, mid-latitude 
erythermal doses in the Southern Hemisphere axe comparable to an equatorxjard 
latitudinal shift of about 7°. In Melbourne, which is somewhat on the equa- 
torward side of the central mid-latitudinal belt, Che suiamer sunburn potential 
under these perturbed conditions would be roughly equivalent to that of pre- 
sent day New Caledonia or along central sections of the Great Barrier Reef. 
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Two additional model calculations have been made concerning the pos- 
sible effects on the 0^ layer due to operation of an SST fleet of aircraft. 

On of these (Run 19) assumes that the SST corridor is centered at an altitude 
of 17km (rather than 20km as in Run 18) at 45 “N, x^hile the other calculation 
assumes a tropical flight corridor at 201on altitude and lO^N. Summer and 
mnter seasonally averaged distributions of total columnar 0^ as functions 
of latitude for these two experiments (Runs 19 and 20) , as well as for the 
unperturbed (Run 17) and the 20km, 45°N SST corridor (Run 18) cases are 
compared in Fig. 12. The calculations show that for the conditions of Run 
19 (corridor at 17km, 45“N) , approximately 6 % of the total unperturbed 0^ 
will be destroyed during an annual cycle. Most of this depletion occurs 
in the Northern Hemisphere where 8% of the total 0^ contained there will 
be lost. In contrast, only about 3% depletion will occur in the Southern 
Hemisphere. This is particularly evident during the Northern Hemisphere 
winter and summer seasons shoTvn in Fig. 12 during v?hich 11% (N. H. winter) 
and 8% (N. H. summer) 0^ depletions are indicated. At the same time, neg- 
ligible 0^ losses are noted in the Southern Hemisphere. Thus, we see, in 
comparison with the results of Run 18, that lowering the altitude at which 
a fleet of SST's would operate reduces its impact on 0^ depletion in the 
Southern Hemisphere, but still must be considered as a factor in the Northern 
Hemisphere . 

The situation is somewhat different, however, when SST's fly in the 
tropics, as simulated in Run 20 (corridor at 20 km, 10°N) . For this case 
approximately 12% of the total global 0^ will be destroyed throughout the 
year (the same amount as obtained from Run 18), but this loss is fairly 



evenly distributed between the two hemispheres. This feature can be seen 
in Fig. 12, especially for the Northern Hemisphere summer during which 
the difference between the unperturbed curve (Run 17) and the Run 20 per- 
turbed curve is fairly uniform from pole to pole. 

The odd nitrogen injection rate assumed in our model calculations (1.8 x 
10 metric tons per year) was chosen to most nearly represent the effect of 
a rather large fleet (about 500 aircraft) of the proposed Boeing 2707 SST, 
which might have existed near the year 2000. However, if, in fact, only 
SST’s of the current Anglo-French Concorde and Russian Tupolev 144 proto- 
types, which are smaller and fly at lower, less .harmful altitudes (^17 km), 
are built, we estimate that it would require a world wide fleer of a few thou- 
sand aircraft to attain an effective inj action rate as large as the one used 
in our model simulation. Clearly, the nature and proliferation of the global 
SST fleet of the future must be carefully analysed if we intend to maintain 
the environment as we Icnox^r it today. 
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Fig. 1 


Fig. 2 
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Fig. 4 


Fig. 5 


Fig. 6 


Fig. 7 


Fig. 8 


Fig. 9 


. Model surface topography ia dekameters drawn from the truncatea 
spectral representation. 

t 

. The atmospheric temperature distribution i' ':ained for season 10, 

Run 12 compared against the observed atmospheric temperature dis- 
tribution for summer and winter according to Newell (1969) . 

. Vertical cross section of the zonal wind (m/s) for season 10, 

Run 12 compared ^rlth the summer-winter observations analyzed by 
Newell (1969). 

12 

. Mean meridional circulation (mass flow in units of 10 g/s) 

redrawn from an analysis by Louis (1974) for the months December- 
February. 

. Mean meridional circulation patterns from the model results 
(Run 17, season 12) for the months December-February. 

, The unperturbed winter-summer NO 2 distribution (in parts per 

billion by volume) according to Hesstvedt (1974) . 

11 

. Ozone number density (10 /cm) cross sections for the winter-summer 
season predicted by the model for season 12, Run 17 compared with 
observations analyzed by Wu (1973) . 

. Columnar 0^ distribution (in Dobson units) in the Northern Hemisphere 
as a function of season and latitude: (a) as determined by Dutsch 
(1971) from observations and (b) as predicted during year 3 (Run 17) 
in our unperturbed model. 

. The winter- summer SST induced perturbation in NO^ (in parts per 
billion by volume) from the two-dimensional model of Prinn, et al 
(1975). 
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T?ig. 10. 


Fig. 11, 


Fig, 12. 


Seasonally averaged distributions of columnar 0^ (in Dobson units) 
as a function of latitude predicted for the unperturbed (dashed 
lines - Run 17) and perturbed (solid lines - Run 18) stratospheres 
Approximate daily mean eryth^^al doses predicted from the model 
Og results shot-m in Fig. 10 but limited to the months of January 
and July. 

Same as Fig. 10 but adding the results of the SST hypothetically 
perturbed conditions of Runs 19 and 20 (see text for details) . 
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